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The orthorhombic boride crystal family XYB14, where X and Y are metal atoms, plays a critical 
role in a unique class of superhard compounds, yet there have been no studies aimed at under- 
standing the origin of the mechanical strength of this compound. We present here the results from 
a comprehensive investigation into the fracture strength of the archetypal AlLiBi4 crystal. First- 
principles, ab initio, methods are used to determine the ideal brittle cleavage strength for several 
high-symmetry orientations. The elastic tensor and the orientation-dependent Young's modulus 
are calculated. From these results the lower bound fracture strength of AlLiBi4 is predicted to be 
between 29 and 31 GPa, which is near the measured hardness reported in the literature. These 
results indicate that the intrinsic strength of AlLiBi4 is limited by the interatomic B-B bonds that 
span between the B layers. 



The development of new superhard materials that can 
operate under extreme conditions is critical for high- 
performance industrial manufacturing and is a subject 
that has recently received great attention. [ll-Q The or- 
thorhombic borides, formulated as XYB14 where X and 
Y are metal atoms, have been of interest to scientists and 
engineers for the past decade due to a report 0, [l] that 
AlMgBi4 prepared by mechanical milling can achieve a 
hardness between 32 and 46 GPa. The reason for the ob- 
served superhardness is not understood. It is suggested 
that in part the strength is due to the so-called "nano- 
composite" microstructure comprised of AlMgBi4 and 
TiB2, although the hardness of each individual phase is 
beheved to be less than the hardness of the composite, m 
There have been many studies examining TiB2, |7H10l| 
but the orthorhombic boride family has received much 
less attention and is therefore the focus of this letter. 

Whereas most hard materials are dense, highly- 
symmetric crystals, the XYB14 structure, shown in 
Fig. [1] is relatively open and has low symmetry (Imma) . 
This crystal structure, which was first reported by 
Matkovich and Economy in 1970, has a unit cell contain- 
ing four formula units of XYB14. [llj The 64-atom unit 
cell can be expressed as X4Y4(Bi2)4Bg to distinguish the 
two B allotropes. The B layers are constructed from B12 
icosahedron that are connected to each other through 
the so-called inter-icosahedra B atoms that are trigo- 
nally bonded to three B12 units within the B layer. Re- 
cent spectroscopy evidence indicates that the B-B bonds 
that span between the B layers, directly connecting icosa- 



hedron, are very strong. 1^ Unlike many metal-boride 



compounds the metal atoms are not covalently bonded 
to the B, but instead the metal atoms ionize and do- 
nate their valence electrons to the covalently bonded B 
network. As a near-superhard material, this crys- 

tal family is unique, which has led us to investigate the 
bonding in the crystal as it affects the crystal's mechan- 
ical strength. 
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FIG. 1. (color online) A simple schematic of the XYB14 crys- 
tal structure. The red (medium grey) spheres are the X-site, 
the blue (dark grey) spheres are the Y-site, and the light grey 
spheres are the B atoms. The short-dashed arrows denoted 
by Ejjjjjj and Emax show the loading directions that yield the 
minimum and maximum Young's modulus for AlLiBi4. The 
six planes selected for examination within the brittle cleavage 
model are shown as long-dashed lines and are labeled accord- 
ing to the nomenclature used in Table H] 



diffraction methods to refine the crystallographic data. 
For some of the compounds, such as AlMgBi4, a rela- 
tively high concentration of vacancies, around 25%, are 
found at the metal atom sites. [l3] Diffraction results 
for other XYB14 compounds, such as AlLiBi4, do not 
find such a large number of vacancies. 18|, |l9[ Wer- 
heit et al. have used Raman spectroscopy to compared 



Following Matkovich and Economy's work,[ll[ Higashi 
and Ito synthesized several XYB14 compounds and used 



the vibrational spectrum of various XYB14 compounds 
and have found that the AlLiBi4 crystal has less inter- 
nal distortions than many other XYB14 structures, in- 
cluding AlMgBi4. 1^ From these results we conclude 
that experimental specimens of AlLiBi4 are likely to have 
fewer point defects than many other XYB14 compounds 
and consequently AlLiBi4 is selected as the archetypical 
structure for study in this letter. 

Previous theoretical studies have focused on the effects 
of chemical substitution on the properties of the XYB14 
crystal family. This is in part because the sys- 
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tern is known to accept a large number of metal species 
including Li, Be, Na, Mg, Al, as well as a variety of 
rare-earth elements, such as Tb, Dy, Ho, Er, Yb, and 
Lu. 0, [11 [11 [2l|-[2l In addition the superhard Ames 
Lab specimen was synthesized by a mechanical alloy- 
ing method, which introduces a wide variety of impurity 
species to the crystal including Ti, Si, Fe, O, and C. In 
these theoretical studies the reported figure of merit for 
hardness is the bulk modulus because its computation 
is relatively simple. However, the bulk modulus alone 
only gives information about the average bond strength 
under an applied volume dilation and docs not give any 
information about the strength of individual bonds in 
the crystal. Understanding hardness requires knowledge 
about the local mechanisms for bond breaking as it re- 
lates to fracture in the crystal. In this letter the frac- 
ture strength of AlLiBi4 will be examined using an ideal 
brittle cleavage model. This approach allows for insight 
regarding the local bonding within the crystal and may 
lead to a strategy for improving the hardness. 

The first-principles, density functional theory method 
used in this study is implemented in the SIESTA soft- 
ware package. [H [2^ The Perdew-Burkc-Ernzerhof gen- 
eralized gradient approximation (GGA) is used for the 
exchange-correlation energy and norm-conserving pseu- 
dopotentials are used in place of the all-electron atomic 
potentials. 27, [51] The wavefunction is represented by 
a set of finite-range numerical atomic orbitals. Each 
atomic basis is extended to include double—^ functions 
plus a shell polarization that is constructed using the 
split-valence scheme. [1^ The cutoff radii used for each 
C function are presented in Ref. [s^. Real space mesh- 
ing is performed to an energy cutoff of 175 Rydberg. 
The Kohn-Sham energies are sampled across the Bril- 
louin zone using al2xl2xl2 Monkhorst-Pack grid. 31 1 
The atomic structural optimization follows the conju- 
gate gradient minimization method and the thresholds 
for the residual forces on atoms and the supercell are 
0.005 eV/A and 0.0005 eV/A^ respectively. The calcu- 
lated lattice parameters for AlLiBi4 are 5.88, 10.39, and 
8.15 A, which agree well with the reported experimental 
values 5.847, 10.354, and 8.143 A. [l| 

The ideal brittle cleavage model used here separates 
the AlLiBi4 crystal into two semi- infinite, rigid atomic 
blocks that are pulled apart to introduce a pair of cleav- 
age surfaces at a predefined atomic plane. This ideal- 
ized approach simultaneously stretches and breaks all the 
bonds at the interface. Although the effect of crack tip 
initiation and propagation cannot be included using this 
method, it allows for the bond strengths localized in the 
crystal to be investigated. Internal atomic relaxations 
and lattice contractions perpendicular to the direction 
of elongation are forbidden. These constraints allow the 
strength of the bonds across the cleavage interface to be 
determined independent of possible near-surface atomic 
reconstructions, which would be present in an experi- 



TABLE I. The numerical results for applying the ideal brittle 
cleavage model to the six AlLiBi4 cleavage planes that are 
shown in Fig. [1] 



O rient cit ion 


Cleavage energy 


Critical lenf 


rth Critical stress 


Gt/Area (J/m^) 


lb (A) 


at (GPa) 


{100}-B 


7.94 


0.51 


57.7 


{100}-M 


7.74 


0.45 


63.5 


{010}-B 


9.16 


0.49 


68.4 


{010}-M 


5.42 


0.56 


35.7 


{001}-B 


8.27 


0.49 


62.0 


{001}-M 


7.51 


0.57 


48.1 



mental specimen. The calculated energy of the cleaved 
crystal, relative to the energy of the perfect crystal, is 
called the decohesive energy, Eb, and is determined as a 
function of the interplaner spacing, x, across the speci- 
fied cleavage planes. The decohesive energy is fit to the 
universal binding energy relation (UBER) developed by 
Rose et al. in Ref. l32l which is expressed more precisely 
in Ref. |33| as. 
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exp 



When all the atomic bonds that span the cleavage in- 
terface are broken, the decohesive energy saturates to 
the cleavage energy, G&. The cleavage stress is the first 
derivative of the decohesive energy with respect to x. 
The critical cleavage stress, ab, is defined as the maxi- 
mum stress, and the corresponding interplanar spacing 
is referred to as the critical length, lb- 

For AlLiBi4 cleavage is considered within the high- 
symmetry {100}, {010}, and {001} families of planes. 
For each crystallographic direction two cleavage planes 
are examined: one that passes through the icosahedron, 
labeled "B," and one that passes between the icosahe- 
dron, labeled "M." These were selected to best represent 
the variation in the bonding for each of the sampled di- 
rections; one of the cleavage planes has many bonds that 
span the interface and the other few bonds. It is intended 
to test interfaces with the highest and lowest fracture en- 
ergies. These planes are identified in Fig. [1] A 128-atom 
supercell is used to guarantee that the calculated deco- 
hesive energies are converged to better than 0.005 J/m^. 
The decohesive energies, UBER fit, and derived stresses 
are plotted in Fig. [2l The computed data matches the 
functional form of the UBER relation very well and the 
resulting critical parameters are listed in Table HI 

The decohesive energy curves and stresses for the 
{100}-B and -M planes are very similar, the critical 
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FIG. 2. (color online) The ideal brittle cleavage results for 
AlLiBi4. The left ordinate axis labels the energies and the 
right the stresses. The top frame gives the results for the 
{100} planes, the middle the {010} planes, and the bottom 
the {001} planes. In each frame the "B" results are red (light 
grey) and the "M" results are blue (dark grey). The decohe- 
sive energy DFT data are solid symbols, the UBER relations 
are solid lines, and the stresses are dashed lines. 



stresses differ by less than 10%. This is not surpris- 
ing considering that the density and geometric arrange- 
ments of the B-B bonds in these planes are nearly equiv- 
alent. Whereas the {100} planes are very similar, the 
{010} planes are considerably different. The {010}-B 
plane passes through the B layer bisecting the icosahe- 
dron, breaking many B-B bonds, but the {010}~M plane 
passes between the B layers and therefore cuts signif- 
icantly fewer bonds. Within the ideal brittle cleavage 
model the calculated critical stress for the M plane is 
48% smaller than that of the B plane. 

According to the Raman spectroscopy results reported 
in Ref. |l2|, the B~B bonds that span between the B lay- 
ers and connect the icosahedron are expected to have a 
greater binding strength than the bonds inside the icosa- 
hedron. This can be examined qualitatively by plotting 
the bonding charge density, as shown in Fig. |31 |34| The 
bonding process results in a buildup of charge in the B-B 
bonds that bridge the B layers. In Fig. |3] charge accumu- 
lation between the B layers is observed both between the 




FIG. 3. (online color) An isosurface showing regions in an 
XYBi4 crystal with positive bonding charge density. The 
atomic sites are color coded following the description in Fig.[l] 



icosahedron and the inter- icosahedra B. The results from 
the ideal cleavage model, presented in Table HI indicate 
that regardless of the anticipated high strength of the 
bonds at this plane, the relatively low number density 
of bonds causes the {010}-M plane to have the lowest 
critical stress of all the planes examined here. It can be 
concluded that the B-B bonds that span between the B 
layers are key for controlling the overall strength of the 
crystal. 

It is the weakest planes that are of primary interest 
because fracture naturally transverses the weakest path 
through a crystal. It is these same planes that also dom- 
inate the elastic response. For a given family of planes, 
the elastic response to a uniaxial load applied normal to 
the planes should be an indicator of the relative cleav- 
age strength, i.e., the Young's modulus should scale with 
the cleavage strength. For the orthorhombic XYBi4-type 
crystal there are nine unique tensor elements that can 
be derived from the linear stress-strain relation and the 
crystal symmetries. The components of the stiffness and 
compliance tensors for AlLiBi4 are calculated and pre- 
sented in Table HIl In Fig.Sl^a), the Young's modulus, i5, 
is represented as a function of crystallographic orienta- 
tion, according to the formula, 
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l-ySw + I2S22 + ^3533 + 21112812 + 2lj^l^si3 + 2I2I3S23 

+ I2I3S44 + ^1^3555 -I- ^1^2*66, 



where are the elastic compliance tensor components 
and h, I2, and ^3 are the direction cosines. The represen- 
tation surface in Fig.HJ^a) is projected on the (100), (010) 
and (001) planes, and the results are shown in Fig. |3^b). 

For the [010], [001], and [100] directions the Young's 
modulus is 293, 404 and 505 GPa. From the ideal brittle 
cleavage model the minimum critical stresses for these 
same directions are 35.7, 48.1, and 57.7 GPa, as listed 
in Table HI Comparing these numbers demonstrates that 
indeed the directional representation of the Young's mod- 
ulus is an accurate predictor of the relative ideal fracture 
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TABLE II. The elastic tensor components for AlLiBi4. 



Stiffness coefficients 


Cll 


C22 


C33 


C44 


C55 


C66 


C12 


Cl3 


C23 


(GPa) 


526 


411 


419 


91.0 


201 


130 


45.7 


83.4 


32.0 


Compliance coefficients 


Sll 


S22 


S33 


S44 


S55 


S66 


S12 


Sl3 


S23 


(xlO-'^mVN) 


1.98 


2.47 


2.47 


10.99 


4.98 


7.69 


-0.19 


-0.38 


-0.15 




FIG. 4. (color online) The representation surface for the 
Young's modulus of AlLiBi4 (given in GPa) as a function 
of crystallographic orientation. Frame (a) shows a 3D plot of 
the representation surface. Frame (b) shows 2D projections 
of the representation surface onto the (100), (010) and (001) 
planes and are printed in green (light grey) , blue (dark grey) , 
and red (medium grey) respectively. 



strength of a particular orientation. From the results 
in Fig. |4] the minimum value of the Young's modulus is 
256.0 GPa, which corresponds to a uniaxial load orien- 
tated (0 = 90°, e = 44.96°), where (f) is the angle of ro- 
tation from the positive x— axis to the positive axis 
in the xy— plane and 6 is the out-of-plane angle of ro- 
tation from the positive z— axis to the xy— plane. This 
direction is drawn in Fig. [T] as a short-dashed arrow that 



is labeled 



Assuming linear proportional scaling. 



the computed results for the high-symmetry orientations 
can be used to predict that the cleavage strength for a 
uniaxial load applied in the Ej^-^jj^^ direction is between 
29 and 31 GPa. This is the predicted lower limit of the 
ideal brittle cleavage strength for AlLiBi4. We submit 
that for a brittle material, such as AlLiBi4, which does 
not undergo extensive plastic deformation near the crack 
tip prior to fracture, the calculated ideal brittle cleav- 
age strength is a reasonable estimation of the fracture 
strength. The physical features of fracture neglected in 
this ideal brittle cleavage model, including crack tip plas- 
ticity, lattice contractions, and atomic reconstructions, 
results in an overestimation of the cleavage energy and 
subsequently the actual critical energy and stress will be 
lower than our calculated results. The experimentally 
measured hardness for AlLiBi4 ranges between 20 and 



29 GPa 3 21 1, which suggests that for this material the 
atomic scale behavior that we have reported here plays 
an important role in determining the actual hardness of 
the material. 

In summary, we have coupled the results from a series 
of ideal brittle cleavage strength calculations to the calcu- 
lated orientation dependent Young's modulus to predict 
the fracture strength of AlLiBi4. While admittedly this 
simplistic model neglects some of the macroscopic fea- 
tures of fracture associated with crack-tips, lattice plas- 
ticity, and interface reconstructions, we believe that here 
we have demonstrated that this still may be an effective 
approach to gauge the strength of brittle materials, such 
as the XYB14 crystal family. In contrast to all of the 
previous theoretical studies of the XYB14 crystal family, 
which have used the bulk modulus as an indicator of the 
bond strength, the approach used here allows for the lo- 
cal bond strength to be investigated on a planc-by-plane 
basis. Unlike the more sophisticated, multi-scale model- 
ing approaches, which have been deployed to study frac- 
ture in p olycrys talline diamond. Si, and other metallic 



systems, [35 



-38| the method used here is relatively sim- 
ple. We believe that our approach can be used to screen 
prospective structures prior to their being investigated 
using a more elaborate theoretical technique. 

The existing picture of bonding in the XYB14 crys- 
tal family is that B forms a covalent network of atoms 
constructed of B12 icosahedron. The B12 are stabilized 
by the electrons donated by the ionized metal atoms, ac- 
cording to the Jemmis mno rules, [s^ Excess charge accu- 
mulates in the inter-icosahedra bonds, both those within 
the B layer and those that span between the layers. Ex- 
perimental results indicate that the inter-icosahedra B-B 
bonds spanning between the layers are stronger than the 
intra-icosahedra bonds, [l^l Here we find that regardless 
the strength of the inter-icosahedra bonds the fracture is 
significantly more likely to proceed between the icosahe- 
dron rather than through, due to the density of bonds at 
the cleavage plane. In fact the {010}-M planes are the 
weakest of those examined in this study, which suggests 
that the hardness of the material may be closely tied to 
the B-B bonding that connects the icosahedra layers. In 
practice, this means the intrinsic strength of this crystal 
family possibly can be changed, either strengthened or 
weakened, by the introduction of a dopant species that 
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directly affects tlicse bonds. 
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